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Abstract 

We compute the chiral logarithmic corr ections to the B% - B d and - B s mixing 
amplitudes in the Standard Model and beyond. We then investigate the impact of 
the inclusion of the lowest-lying scalar heavy-light states to the decay constants and 
bag-parameters and show that this does not modify the pion chiral logarithms, 
but it does produce corrections which are competitive in size with the K- and ij- 
meson chiral logarithms. This conclusion is highly relevant to the lattice studies 
since the pion chiral logarithms represent the most important effect in guiding the 
chiral extrapolations of the lattice data for these quantities. It is also important to 
stress that the pion chiral logarithmic corrections are useful in guiding those chiral 
extrapolations as long as m n <C As, where As stands for the mass difference between 
the heavy-light mesons belonging to ^ + and ^ doublets. 



1 Unite mixte de Recherche du CNRS - UMR 8627. 



1 Introduction 



The oscillations in the B% s — B ds systems are mediated by the flavor changing neutral 
currents which are forbidden at tree level of the Standard Model (SM) and therefore their 
detection gives access to the particle content in the corresponding loop diagrams. First 
experimental measurement of a large value of Am^ indicated that the top quark mass 
was very heavy [I], which was confirmed almost a decade later in the direct measurements, 
m t = 172.5 ± 1.3 ± 1.9 GeV through the pp-collisions [2]. Nowadays, the accurately mea- 
sured AniB d = 0.509(5) (3) ps" 1 [3J, and Am Bs = 17.31 (fy) (7) ps -1 [3], are used to constrain 
the shape of the Cabbibo-Kobayashi-Maskawa (CKM) unitarity triangle and thereby de- 
termine the amount of the CP- violation in the SM [5]. This goal is somewhat hampered 
by the theoretical uncertainties in computing the values for the two decay constants, /s sd , 
and the corresponding "bag" parameters, Bs 3d - These quantities can, in principle, be com- 
puted on the lattice. 1 However, a major obstacle in the current lattice studies is that the 
c?-quark cannot be reached directly but through an extrapolation of the results obtained 
by working with larger light quark masses down to the physical c?-quark mass. This ex- 
trapolation induces systematic uncertainties which are hard to control as the spontaneous 
chiral symmetry breaking effects are expected to become increasingly pronounced as one 
lowers the light quark mass [7|. Heavy meson chiral perturbation theory (HMChPT) allows 
us to gain some control over these uncertainties because it predicts the chiral behavior of 
the hadronic quantities relevant to the heavy-light quark phenomenology which then can 
be implemented to guide the extrapolation of the lattice results. HMChPT combines the 
heavy quark effective theory (HQET) with the common pattern of spontaneous breaking 
of the chiral symmetry, SU(3) L ® SU(3) R -> SU(3) V 0. 

Like in the standard ChPT, in HMChPT one computes the chiral logarithmic corrections 
(the so-called non-analytic terms) which are expected to be relevant to the very low energy 
region, i.e., m q <C Aqcd- While this condition is satisfied for u- and <i-quarks, the situation 
with the s-quark is still unclear [9]. Also ambiguous is the size of the chiral symmetry 
breaking scale, A x . Some authors consider it to be around A7rf n ~ 1 GeV [TU], while others 
prefer identifying it with the mass of the first vector resonance, m p = 0.77 GeV (see e.g. 
ref. [HJ), and sometimes even lower [12] . In the heavy-light quark systems the situation 
becomes more complicated because the first orbital excitations (j[ = l/2 + ) are not far 
away from the lowest lying states (j[ = l/2~). The recent experimental evidence for the 
scalar Dq s and axial D ls mesons indicate that this splitting is only Ag s e mo^ - m D a — 
m,£) ls — mo* = 350 MeV [IB], and somewhat larger for the non-strange states A$ ud = 
430(30) MeV [2]. 2 This and the result of the lattice QCD study in the static heavy 
quark limit [T7j suggest that the size of this mass difference remains as such in the 6-quark 
sector as well. One immediately observes that both As B and A$ u d are smaller than A x , m v , 
and even mjc, which requires revisiting the predictions based on HMChPT and reassessing 
their range of validity. In this paper we investigate this issue on the specific examples of 

1 Recent reviews on the current status of the lattice QCD computations of — ~lf q mixing amplitudes 
can be found in ref. [6] . 

2 We note, in passing, that the experimentally established fact, As s < Ag u d , is not yet understood [15] 
although a recent lattice study with the domain wall quarks indicates a qualitative agreement with exper- 
iment |16j . 
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the decay constants /r, s and the bag parameters which enter the investigation of the SM 



and supersymmetric (SUSY) effects in the B® — B° d and B® — T3" s mixing amplitudes 

2 Bases of AB = 2 operators and B-parameters 



=r0 



The SUSY contributions to the B q — B q mixing amplitude, where q stands for either d- or 
s-quark, are usually discussed in the so called SUSY basis of AB = 2 operators 



Oi = 




- 75 )g^V(i-75)g J 


o 2 = 


b\l - 


75 y^(l- 75 y, 


o 3 = 


6*(1 - 


7 5 V^(i-75)g\ 


4 = 


b\l- 


75 )g l ^(l + 75 )g J , 


o 5 = 


b l {l- 


J 5 )q j V(l + j 5 )q\ 



(1) 



where i and j are the color indices. Although the operators in the above bases are written 
with both parity even and parity odd parts, only the parity even ones survive in the matrix 
elements. In SM, only 0± (left-left) operator is relevant in describing the B q — B mixing 
amplitude. The matrix elements of the above operators are conventionally parameterized 
in terms of bag-parameters, Bi_ 5 , as a measure of the discrepancy with respect to the 
estimate obtained by using the vacuum saturation approximation (VSA), 

mo^iB*) 

Bi- 5 {v) , (2) 



(ajlOi-sOWvaA 



where v is the renormalisation scale of the logarithmically divergent operators, 0«, at which 
the separation between the long-distance (matrix elements) and short-distance (Wilson 
coefficients) physics is made. We remind the reader that 



(SJIOilfljW = 2(1 + -)<B 9 VmIO> <0|^|B 9 > 



(s;io2|s;>vsa = -2(i-i) |(o|p|5°)| 2 



m0 3 \B») YSA = \l- 2 -) \(0\P\B°\r , (3) 
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(B° q \O A \B q ) VSA = \(B° q \A,\0) (0\A»\B° q )+2\{0\P\B q >)\ 2 , 
(B° q \0 5 \B° q ) VSA = (B° q \A,\0) (0\A»\B° q } + 2 -\(0\P\B° q }\\ 

with = brf^^q and P = bj^q being the axial current and the pseudoscalar density, 
respectively. In HQET, in which we will be working from now on, the field b is replaced by 
the static one, h\ which satisfies ftJ 7o = h/. This equation and the fact that the amplitude 
is invariant under the Fierz transformation in Dirac indices, eliminate the operator O3 from 
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further discussion, i.e., (B q \Os + O2 + \Oi\B q ) = 0, where the tilde is used to stress that 
the operators are now being considered in the static limit of HQET (\v\ = 0). Furthermore, 
in the same limit 

hm -£= = hm -== = (0 4, .B£(i>))hqet = if q , (4) 

where f q is the decay constant of the static l/2~ heavy-light meson, and the HQET states 
are normalized as (B^(v)\B^(v')) = 5(v — v'), so that we finally have 

= lm 2 B lq (u), 

(B q \O 2 (u)\B° q ) = - Z -f q {vfB 2q {v), (5) 
(B° q \d 4 (u)\B° q ) = 7 -f q (v) 2 B 4q (v), 
(B° q \d 5 (u)\B° q } = b -f q {vYB bq {v). 

One of the reasons why lattice QCD is the best currently available method for comput- 
ing these matrix elements is the fact that it enables a control over the //-dependence by 
verifying the corresponding renormalisation group equations, which is essential for the can- 
cellation against the ^-dependence in the corresponding perturbatively computed Wilson 
coefficients [20J. From now on we will assume that the UV divergences are being taken care 
of and the scale v will be implicit. 



3 Chiral logarithmic corrections 

In this section we use HMChPT to describe the low energy behavior of the matrix ele- 
ments ([5]). Before entering the details, we notice that the operators O4 and O5 differ only 
in the color indices, i.e., by a gluon exchange, which is a local effect that cannot influence 
the long distance behavior described by ChPT. In other words, from the point of view of 
ChPT, the entire difference of the chiral behavior of the bag parameters B iq and B 5q is 
encoded in the local counter-terms, whereas their chiral logarithmic behavior is the same. 
Similar observation has been made for the operators entering the SUSY analysis of the 
K°-K° mixing amplitude, as well as for the electromagnetic penguin operators in K — ► tttt 
decay [21]. Thus, in the static heavy quark limit {ttiq — > 00), we are left with the first three 
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operators in eq. (J5j) which, in their bosonised version, can be written as [22] 
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^/3 lx Tr (£ir\ 7/ ,(l - 7s )X Tr {^H%^{l-^)X 



x 



c.t. 



]T/? 2X Tr (£ir\(l- 75 )X Tr (£H%(1 - 7s )X 



+ c.t. 



]T/3 4 xTr (£iT\(l- 75 )X Tr (^) ? (1 + 75 )X 



A' 



+/3 4X Tr - 7s )X Tr (^H^l + ^X 



+ c.t. 



(6) 



where X G {1, 7 5, 7i ,, J u j5, &v P }- 4 As before the index "g" denotes the light quark flavor, 
and "c.t." stands for the local counter-terms. To relate /3j's to the bag parameters in eq. 
we should recall that the field H q is built up from the pseudoscalar (P) and the vector (P*) 
meson fields as 



! + {> 



[pQ*( v )^- P Q( v ) l5 ] 



pQ*( v )Y-P Q (v) lt 



1 - 



(7) 



thus exhibiting the heavy quark spin symmetry for the lowest lying j[ = 1/2 states. After 
evaluating the traces in eq. (jSJ) we obtain 



By 



2f 



-A 



Bo 



12 

5f 



12 
7? 



1 • 



12 

5? 



4 • 



(8) 



where 



Si 



A + /3i 7B - 4(/3i > + /3i 7l/76 ) - 12/3^, 

— $2 — ^275 + + A^tb, 



(9) 



3 In the first version of the present paper only the factorisable pieces in the bosonised forms of the 
operators 02,4 were considered. We thank David Lin for pointing out to us the presence of the corresponding 
nonfactorisable pieces, recently considered in rcf. |22| . which are being properly accounted for in this version 
of our paper. 

4 Contraction of Lorentz indices and HQET parity conservation requires the same X to appear in both 
traces of a summation term. Any insertions of jj can be absorbed via jjH = H, while any nonfactorisable 
contribution with a single trace over Dirac matrices can be reduced to this form by using the 4x4 matrix 
identity 

4Tr(AB) = Tr(A)Tr(B) + Tr( 75 A)Tr( 75 B) + Tr(A 7AI )Tr( 7 ' J S) 
+Tr(^ 7M75 )Tr( 757 ^S) + 1 / '2Tr (Aa ^)Tr (a^ B). 
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We will use the well known form of the HMChPT lagrangian [8] 

£ = flight + Cl- + £ct. j 

2 

flight = ^tr [(9 M St) + Et x + X + E ] , 

o 

Ci- = iTr [H b vD ba H a ] + gTr [H blfll5 Al~H a ] , 

£ ct = hTr [H a H b (£MZ + e^e^J + hTr \H a H a (£M£ + ?M?) bb ] , (10) 

where H a (v) = JqHKv)^, g is the coupling of the pseudo-Goldstone boson to the pair of 
heavy-light mesons, and 

DlH b = d»H a - H b h&%t + CdAa , Af = - £8^% , 

£ = S = exp fay) > .M = diag(m„,m d ,m s ) , 



V 



2" ^75^ 
K~ 



7T+ AT+ \ 

* / 



11^ 



with / m 130 MeV, and x = 2B A4. In the above formulae H a refers to either or H®, 
defined in eq. ((?]). Note also that we distinguish between the trace over Dirac ("Tr") and 
flavor ("tr") indices. In the chiral power counting the lagrangian £n g ht in eq. (ITU!) is of 0{p 2 ) 
while the rest of £ is of 0{p 1 ). To get the chiral logarithmic corrections to -Bjg-parameters, 
we should subtract twice the chiral corrections to the decay constant f q from the chiral 
corrections to the four-quark operators ([6]). The former is obtained from the study of the 
bosonised left-handed weak current 



(V-AY 



i a 



{Tr [(ftf ), 7 "(1 - 7s)] (1 + xfaM) + x x Tr [^MH) q Y(l - 7s)] } , (12) 



where a is the tree level decay constant in the chiral expansion, and x 12 are the counter-term 
coefficients. Together with the strong coupling g, these parameters are not predicted within 
HMChPT. Instead, they are expected to be fixed by matching the HMChPT expressions 
with the results of lattice QCD for a given quantity (see reviews in ref. [23J ) . The notation 
used above is the same as in ref. [21]. The chiral logarithmic corrections to the decay 
constant come from the diagrams shown in fig. [1] 

1 /3 „ m 2 1 „ m?, 1 „ rr 2 

1 



fd 



a 



(47T/) 5 



3 2 . m- 1 o . m K 
-m n log — + -m K log ■ 

4 a 2 



1 2 m 2 \ 
12^ l0g ^J 



+x 1 (/i)md + x 2 (/i)(m u + m d + m s ) + -5Z d 



fs 



a 



(4tt/) ! 



2 i m /C 1 2, 

m ^ log ^ + 3^ log 7^J 



(13) 



+xi(/i)m s + x 2 (fi)(m u + m d + m s ) + 
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\ 

\ 
I 



o- i- o- o- 

Fi gure 1: The diagrams which give non-vanishing chiral logarithmic corrections to the pseudoscalar 
heavy-light meson decay constant. The double line indicates the heavy-light meson and the dashed one 
the pseudo-Goldstone boson propagator. The square stands for the weak current vertex. The full dot is 
proportional to the coupling g. 



i \ 
i \ 
i \ 
i \ 



0" I" I" o- o- o- 

Figure 2: The diagrams relevant to the chiral corrections to the SM bag parameter B\ q . In the text 
we refer to the left one as "sunset", and to the right one as "tadpole". Only the tadpole diagram gives a 
non-vanishing contribution to the bag parameters i?2.4<j- 



where it should be stressed that we work in the exact isospin limit (m u = m d ) so that 
the index d means either u- or c/-quark. Only explicit in the above expressions is the term 
arising from the tadpole diagram (right in fig. [I]), whereas Z^ s , the heavy meson field 
renormalization factors, come from the self energy diagram (left in fig. [1]) and they read 



Zd = l -Wff l0g ^ + mK bg ^ + 6 ^ bg ^ ) 



+ki(n)m d + k 2 {fi){m u + m d + m s ), (14) 
3g 2 ( m? 2 m 2 \ 

{Ait}) 2 V ^ 3 fi 2 J 

In both eqs. (fT3l) and (1T41) the \x dependence in the logarithm cancels against the one in the 
local counter-terms. 

With these ingredients in hands it is now easy to deduce that the only diagrams which 
contribute to the SM bag parameter, _B 1(? , are the two shown in fig. [2j They arise from the 
two terms in O x = A^P;) q ^P^) q + (£P) q (£P) q ) and yield 

"tadpole": - 4^^^ ^(ty 2 m 2 log ^ , (15) 
respectively, where f are the SU(3) generators and rrii masses of the pseudo-Goldstone 
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bosons. The SM bag parameters now read 



D 



id 



Tree 



3s 2 



(4tt/) s 



-77V log + g^log , 



+b 1 {n)m d + b[(fi)(m u + m d + m s 



(16) 



Tree 



3s 2 2 



m 



(4tt/)2 3 



m 2 \og—l + 6i(/x)m s + b[(n)(m u + m d + m s ) 



where we also wrote the counter-term contributions and, for short, we wrote Bj* ce = 
3/3i/2a 2 . The above results agree with the ones presented in refs. [251 [26] . in which the 
pion loop contribution was left out, and with the ones recently presented in ref. [27] . 
As for the bag parameters B 2q and B^ q we obtain 



2Ad 



D Tree 
^2,4 



1 + 



Sg 2 Y T 1 
(4tt/) 2 



ml 



-m„ log 



1 2 

-m log — 



5 



2,4s 



B. 



Tree 
2,4 



2 " ~ /i 

+b 2 ^{p)m d + b' 2A (n)(m u + m d + m 

1 + o 7^ s\9 m v lo S^ + b2AN m s + b 2i {n){m u + m d + m s 
6 {A7rj) z 1 fi z 



(17) 



where B 2 Tree = 12&/4a 2 , £? 4 Tree = 12&/7a 2 Y = @Z t J% A ), with = (3 2lv +(3 2 ^ 5 +A(3 2(7vp , 

and PI = P± lv ~Pa 1v15 + Pa 1u —Pi luTo ■ We checked that our results agree with those presented 
in ref. ^27\ where also the partially quenched theory has been considered. In our paper we 
refer only to the full (unquenched) theory. 



4 Impact of the l/2 + -mesons 

In this section we examine the impact of the heavy-light mesons belonging to the l/2 + 
doublet when propagating in the loops onto the chiral logarithmic corrections derived in 
the previous section. We first extend the lagrangian by adding to eq. (fTOl) the following 
terms [Hj: 

£1+ = -Tr [S b (iv-D ba + A s )S a ) + gTr [S b ~f^ 5 A% a S a ] , 

£ mix = hTr [StftftAfcHa] +h.c, (18) 
C ct = k 1 Tr [S a S b (£M£ + tfMZ^J + h?r [S a S a (£M£ + ?M£%] 

+ k[Tr [H a S b (£M£ + £W) J + k' 2 Tr \H a S a + £W)J + h - c - > 

where the fields of the scalar (Pq) and the axial (Pf ) mesons are organised in a superfield 
S q (v) = [P^(v) W -P (v)] q , = 7 o5j(«)7b. (19) 
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g is the coupling of the P-wave Goldstone boson to the pair of l/2 + heavy-light mesons, 
and h is the coupling of the S'-wave Goldstone boson to the heavy-light mesons, one of 
which belongs to 1/2" and the other to l/2 + doublet. Before including the l/2 + doublet 
we were free to set A = because all the chiral loop divergences are cancelled by 0(m q ) 
counter-terms in the static heavy quark limit. Once the l/2 + doublet is included, the mass 
difference between the l/2 + and l/2~ states (A5 ~ 400 MeV) must be included in the 
lagrangian, but since it does not vanish in the chiral nor in the heavy quark limit it is 
of O(p ) in the chiral power counting (see also ref. [28] where, in addition to the static 
heavy quark limit, the chiral power counting is discussed also when the 1/mQ-corrections 
are included). 



4.1 Decay constants 

Beside the lagrangian, the l/2 + mesons contribution should also be added to the left vector 
current (1121). which now reads 



(V - AX 



icy 

yTr[(^)^(l-75 



■Tr [(6S% 7 M (1 " 75 



+ 



1 a 

T 

ia 



XiTr [(£MH) q >f(l 



7s 



ia 



ia 



+ _x 2 Tr[(^)^(l-75)]tr^ + 



■x 2 Tr[(^) 9 7 At (l-7 5 )] trJW,(20) 



where a + is the coupling of one of the l/2 + mesons to the weak left current, and xi i2 are the 
coefficients of two new counter-terms. From the recent lattice results reported in ref. [29J 
we extract, a + /a = 1.1(2). In other words, at least in the static heavy quark mass limit, 
the weak current coupling of l/2 + mesons is not suppressed with respect to the 1/2 - ones. 
Since we focus on the pseudoscalar meson decay constant, it should be clear that only the 
scalar meson from the l/2 + doublet can propagate in the loop. The diagrams that give 
non- vanishing contributions are shown in fig. [3] and the corresponding expressions now read 



1 + 



(4tt/)= 



>o dx 



J^ml A s ) + J 2 (ml A, 



d 



+A s —(.J 1 (ml A s ) + J 2 (ml A s )) 



a 



b qa b aq 

2(4^7? 



d 



3g 2 lim — [xJiK, x)] - h(mf) - h 2 J^mf, A s ) + J 2 (mf, A s 



*o dx 



( JxK 2 , As) + J2K 2 , A s )) ) - 2h^- (h(m 2 ) + J 2 (m 2 , A s )) 



(21) 



where the summation over "2" is implicit, and we omit the counter-term contributions to 
make the expressions simpler. We stress that in all our formulae, the terms of 0(p 3 ) and 
higher in the chiral power counting are neglected. The integrals and Ji j2 are the same 
as the ones used in ref. [21] and can be found in the appendix of that paper. The terms 
proportional to h 2 comes from the inclusion of the left diagram shown in fig. El while the 
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o- o+ o- o- 0+ 

Figure 3: In addition to the diagrams shown in fig. QJ these two diagrams contribute the loop corrections 
to the pseudoscalar meson decay constant after the 1/2+ mesons are included in HMChPT. The full dot 
in these graphs is proportional to the coupling h. 



last term in the decay constant (the one proportional to h) comes from the right graph in 
fig. [31 Obviously such terms were absent before including the scalar mesons. Notice also 
that 



d 771 ~ 

lim — [xJi(m%, x)] = — m 2 log — I 



*o dx 



(22) 



When As > rrii, which in our case is true for the pion mass, we can expand around m 2 = 
and obtain 



-A, 



A s ) + J 2 (ml, A s ) ■ 
Ii{ml)+I 2 (ml,Ai 

[j l {ml,A s ) + J 2 (m',A< 



2A|(1 - log 



4A|, 
^ > 



ml(l + ]og— £) + ... 



/i 2 



2mt + . 



(23) 



where the dots stand for higher powers in m%. In other words, the presence of the nearby 
l/2 + state does not affect the pionic logarithmic behavior of the decay constant. It does, 
however, affect the kaon and 77-meson loops because those states are heavier than As 
(rriir < As < rriK < rn v ) and the coefficients of their logarithms, although still predictions 
of this approach, cease to be numerically relevant because those logarithms are competitive 
in size with the terms proportional to A| log(4A|//i 2 ), as indicated in eq. (1231) . Stated 
equivalently, the relevant chiral logarithmic corrections are those coming from the SU (2)^® 
SU(2)r — > SU(2)v theory, and the pseudoscalar decay constant reads 




(24) 



where c/(/x) stands for the combination of the counter-term coefficients considered in the 
previous section. 5 At this point we also note that we checked that the chiral logarithms 
in the scalar heavy- light meson decay constant, which has recently been computed on the 



5 More specifically, 2i?oc/(/x) 



ih 2 



[3 + log(4A|/ M 2 )] 



:[l + log(4A|/ M 2 )] = -k!(p) 



4(4tt/) 2 L bv b/M '' 1 2a(47r/) 2L ~ ' '> 2' 

-k[(ii) + k2([i) + k' 2 (fJ,) + >ci{fJ.) + 2x2(11) , where we use the Gell-Mann-Oakes-Renner formula, m 2 = 2Bomd- 
The exact isospin symmetry (m u = rrid) is assumed throughout this work. 
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lattice in ref. [29], are the same as for the pseudoscalar meson, with the coupling g being 
replaced by g, i.e., 



J a 



a 



1 + 3? 2 3 2 



■in 



-m^log— f + ci{fi)m v 



(25) 



Since g^/g 2 ~ 1/9 [30], the deviation from the linear behavior in m\ is less pronounced 
for /+ than it is for f q . Finally, it should be emphasized that the counter-term coefficients 

relevant to the SU(2)y-iheory, namely c^(ju) in eqs. (12411251) . are not the same as those in 
SU{3) V . 

4.2 Bag parameters 

In this subsection we show that the situation with the bag parameters is similar to the 
one with decay constant, namely the pion loop chiral logarithmis remain unchanged when 
the nearby scalar meson is included in HMChPT. To that end, besides eq. ( fl8|) . we should 
include the contributions of l/2 + -mesons to the operators Generically the operators 
0\ 2 4 now become 



Oi 



x 



$>ixTr 
+/3' 1X |Tr 



7 „(1 - To)x\ Tr UtH*) 7 M (1 - 7s)* 



#T 7ft (l-75)X 



Tr 



+ h.c. 



Tr 



7 M (1-75)X 



(26) 



where f3' lx are the couplings of the operator 0\ to both 1/2 and l/2 + mesons, while f}" x 
come from the coupling to the l/2 + mesons only. Similarly, the operators 02,4 now read: 

2 = 5>2*Tr (1- T5 )X Tr ({H^ (1 - 76 )X 

+/3 2X |Tr (pf)^l- l5 )X Tr (^) ? (1 - 7 s)*J + h.c 
+/3 2 ' x Tr (1 - 75)^ Tr [fes*) (1 - 7 8 )*| , (27) 
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o- i-,o+ 0" 0" 0" 0" o- 



i-,o+ i-,o+ 0" 0" 



o- o+ o- o- 0+ o- 

Figure 4: All diagrams which enter in the calculation of the chiral corrections to the operators (Oi^a)- 
Note that we indicate when either the vector or scalar meson can propagate in the loop by 1 _ ,0 + , in 
contrast to the cases when only the scalar (0 + ) meson can couple. 



Oa 



x 

+&x {Tr 
+Kx {Tr 



£iT (l- l5 )X Tr 
(1- 76 )X Tr 

\ (l-Tfe)X 
J q 

(1- 75 )X 



?H Q ) (l + l5 )X 



Tr 
Tr 



(1+7*)* 



+ h.c. 
+ h.c. 





1 (1-75)* 

q 


Tr 




) (1+75)* 




I (1-75)^ 
q 


Tr 




) (1+7 5 )X 
/ q 



(28 



After evaluating the traces in eqs. ( 1261) . and keeping in mind that the external states are 
the pseudoscalar mesons, we have 



Ox = 4(3, 
+0[ 



ep ( 



+ Up up 



epn u p 



+4/?n^) (^p ( 



y ' 

9 



(29) 



where 0, have forms analogous to the ones written in eq. (jSJ). In addition, in eqs. (126f27f28p . 
the fields and S® are defined in a way similar to eq. (J7J). The corresponding tree and 
the 1-loop chiral diagrams are shown in fig. HJ Since the couplings of the four-quark oper- 
ators to the scalar meson are proportional to 0(" and of the pseudoscalar decay constant 
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to a + , the cancellation between the chiral loop corrections in the operators Oi and in the 
decay constant is not automatic. For that reason, instead of writing the chiral logarithmic 
corrections to the bag-parameter, we will write them for the full operator, namely 



-A 



j-i y-it 



-6g 2 lim -^-[xJi(m 2 , x)] + 2I 1 (m 2 ] 
x^o ax 



d 



+2h 2 (J 1 (m 2 , A s ) + Mm*, As) + A s —- (Ji(m 2 , A 5 ) + J 2 (m 2 , A 5 )) 

ai\ s 



+4/4 (h{mt)+h{mlA s )) 
Pi 



h qq h qq 



6g 2 lim ^-[xJi(m 2 , x)] + 2I 1 (m 2 ) 



+4/4 (h{m 2 ) + J 2 (m 2 , As)) + h 2 §- £ J fc (m 2 , sA s )] ) 

Pi Pi fc=l,2;s=±l - 1 J 



(30) 



To keep the above expression simpler we do not write the counter-terms since their structure 
remains the same as before, i.e., a constant times m q and an another constant times tr.M. 
The similar formulae for B 2t 4 q f 2 are lengthy and we will not write them explicitely. For 
the point that we want to make in this section it is enough to consider eq. ( 13"U|) because 
in the expressions for B 2 ^ q f 2 the loop functions / 1>2 and J 12 occur in the same form as in 



eq. (130|) . In the evaluation of the sunset diagrams we used the standard simplification-trick 



(vp - A)(vp -A') A - A' \vp - A vp — A' 



(31) 



We now turn to the case C A5 and study the behavior of eq. (j30~l) around m 2 — > 0. 
In addition to the limits discussed in eq. (123j) . when dealing with the integrals in the last 
line of eq. ( 1301) we shall proceed similarly to what has been done in ref. [31], namely we 
expand the integrand in E n /As and write 



^2 Jk(ml,sA< 



-2(4vr) 2 ^ x ifi e 



c/ 4 £ p 



p^p u 



2(4tt s 



A 



—v IJl v v 



(2tt) 4 ^ (p 2 - ml)[A 2 - (vp) 2 } 
d 4 ~ e p 



p^p v 



(2ir)^p 2 



mi 



0(1/A 2 ) 



2A| 



|(4^)/iK) 



2A 2 10g ^ 



(32) 



where the ellipses stand for the terms of higher order in m 2 /A|. Note, however, that 
even the leading term is already of higher order in the chiral expansion and thus the terms 
proportional to h in eq. fl3Tj|) do not affect the leading chiral logarithmic corrections. 

On the basis of the above discussion and eqs. (I22|I23I) we see that after expanding eq. (IHUl) 
around m 2 = 0, the leading chiral logarithms arising from the pion loops remain unchanged 
even when the coupling to the scalar meson is included in the loops. On the other hand, as 
discussed in the previous subsection, the logarithms arising from the kaon and the ?7-meson 
are competitive in size with those arising from the coupling to the heavy-light scalar meson, 
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which is the consequence of the smallness of As- Therefore, like for the decay constants, 
the relevant chiral expansion is the one derived in the SU (2) l <g> SU (2) r — > SU (2)y theory, 
i.e., 



B lqfq 
B 2Aqfq 



1 - 



m 



(4tt/)s 



1 - 



3# 2 (3-Y) + 3±l 
2(4^7? 

or by using eq. (1241) . for the bag parameters we obtain 



m 



m n log ~f + co 2 Av)m m 



(33) 




(34) 




3g 2 Y T 1 
2(4vr/)2 



m 



^!og-f + CB 2 , 4 (/i)m T 
A 4 



(35) 



which coincide with the pion loop contributions shown in eqs. (1161) and ([IT]) , as they should. 



5 Relevance to the analyses of the lattice QCD data 

It should be stressed that the consequence of the discussion in the previous section is mainly 
important to the phenomenological approaches in which the sizable kaon and ?7-meson loga- 
rithmic corrections are taken as predictions, whereas the counter-term coefficients are fixed 
by matching to large N c expansion or some other model. We showed that the contributions 
of the nearby heavy-light scalar states are competitive in size and thus they cannot be 
ignored nor separated from the discussion of the kaon and/or ?7-meson loops. 

In the extrapolation of the lattice data, instead, this is not a problem because the kaon 
and the ?7-meson loops essentially do not alter the quark mass dependence, whereas the 
important nonlinearity comes from the pion chiral loops. As an illustration, in fig. [5] we 
plot the typical chiral logarithm, — m\ log(m 2 / n 2 ) , as a function of r = rrid/m s which appear 
in the Gell-Mann-Oakes-Renner formulae, 

2 2 r +l2 r+2 

m n = 2B m s r , m K = 2B m s ^— , m v = 2B m s —^— , (36) 

with 2B m s = 2m 2 K — = 0.468 GeV 2 . Thus the fact that the nearby scalar heavy- 
light mesons do not spoil the pion logarithmic corrections to the decay constants and the 
bag-parameters is most welcome from the lattice practitioners' point of view, because the 
formulae derived in HMChPT can still (and should) be used to guide the chiral extrapola- 
tions of the lattice results, albeit for the pion masses lighter than As- 
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Figure 5: Typical chiral logarithmic contributions — mf log(mf //i 2 ) are shown for pion, kaon and n as a 
function of r = m,d/m s , with m s fixed to its physical value, and jj, = 1 GeV. 



6 Conclusion 

In this paper we revisited the computation of the B® — ~B° q mixing amplitudes in the frame- 
work of HMChPT. Besides the SM bag parameter, we also provided the expressions for the 
chiral logarithmic correction to the SUSY bag parameters. More importantly, we study the 
impact of the near scalar mesons to the predictions derived in HMChPT in which these 
contributions were previously ignored. We showed that while the corrections due to the 
nearness of the scalar mesons are competitive in size with the kaon and rj meson loop cor- 
rections, they do not alter the pion chiral logarithms. In other words the valid (pertinent) 
ChPT expressions for the quantities discussed in this paper are those involving pions only. 
This is of major importance for the chiral extrapolations of the results obtained from the 
QCD simulations on the lattice, because precisely the pion chiral logarithms provide the 
important guidance in those extrapolations as long as <C A5. The corresponding useful 
formulas are given in eqs. (I24|34ll35p . As a side-result we verified that the chiral logarithmic 
corrections to the scalar meson decay constant are the same as to the pseudoscalar one, 
modulo replacement g ^>g (c.f. eq. fl^o^) ). 
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